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We discuss modifications of hadron properties in the heat bath extracted from Eu- 
clidian correlation functions and spectral functions reconstructed with the Maxi- 
mum Entropy Method. To investigate the cut-off dependencies we perform simu- 
lation on various lattice sizes and present the result of an analytic calculation of 
the mesonic spectral functions in the infinite temperature limit. 



1. Introduction 

Oh. 

The change of hadron properties in the heat bath is of fundamental interest 
in heavy ion collisions. Especially non-perturbative in-medium effects in the 
vector channel, like the broadening of the resonance peak or a temperature 
, dependence of the mass, could explain the puzzle of the experimentally 

observed low-mass enhancement of the dilepton spectrum. These properties 
of the hadrons are encoded in Euclidian time correlation functions (CF) 
Gh (t, T) which contain in principle all the spectral information in a given 
quantum number channel H . They can be calculated at a discrete set of 
times t € [0, l/aN T ] in the framework of lattice regularized QCD. A more 
direct access to these quantities is provided by the spectral function (SF) 
CTif(w,T) which is related to the correlation function through the integral 
equation 

oo 

Gh(t,T) = J d^(j H ( X )jl(0)) = J ^T) CK ^ { ~iItP ] ^ (1) 

o 

with the renormalized current J^h{t,x) = Zh'4'(t,x)Th4 , ( t t : ^)- The direct 
inversion of this equation to get the SF is an ill posed problem, but became 
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feasible by the Maximum Entropy Method (MEM) 1 . This method allows to 
obtain the most probable spectral function without a priori assumptions on 
the spectral shape. It was successfully applied at zero temperature but its 
application to non-zero temperature due to the limited temporal extension 
of the system remains a challenge. 



2. The Free Spectral Functions 



Up to now MEM has only been applied to data calculated within the Wil- 
son discretization scheme for the fermion action. The corresponding SF 
show a broad bump at high ui which was interpreted 2 as reflecting a bound 
state of a light quark and a heavy Wilson doubler. To investigate this 
explanation in more detail and to study cut-off effects we calculated the 
free spectral function in the chiral and thermodynamic limit. Performing 
a variable transformation 3 and the notation from 4 the free lattice CF for 
Wilson fermions can be written as 



Gfree,lat 
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(2) 



with au{k) = bn + cjj(fc) as channel depending quantities. Therefore all 
cut-off dependencies are encoded in cr^ which is shown for various chan- 
nels in fig.l. They agree with each other at energies below au < 1.5 
which is a remanent of the chiral symmetry, but the explicitly chiral sym- 
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Figure 1. The free spectral function with Wilson-fermions in various quantum number 
channels divided by continuum SF u c = 3/87r _2 w 2 tanh(a;/4T) and with N T = 24. 
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Figure 2. Temperature dependence of the screening masses in various channels (left) 
and the vanishing of the pion pole mass with the quark mass below T c (right). 



metry breaking Wilson term in the fermion action leads to quite differ- 
ent SFs of the chirally symmetric channels (pseudo-scalar(PS)/scalar(S) 
and vector(V)/axial-vector(AV)) at higher energies. The two cusps at 
Loa w 2.2,3.2 and the cut off at uja w 3.9 reflect the three corners of the 
Brillouin zone which restrict the momentum phase space on the lattice. 
Therefore the naive picture that the fermion doublers appear as peaks in 
the SF with masses of 0(1/ a) is perhaps too simple. However, the effect of 
the lattice cut-off on the low energy regime is small. 



3. Temperature dependence of Correlation Functions and 
Spectral Functions 

In the interacting case we have performed quenched simulations with the 
standard plaquette gauge action and non-perturbatively improved clover 
fermions on isotropic lattices with various sizes N% x iV r to investigate 
cut-off and finite volume effects. We simulated up to five different quark 
masses below and in the chiral limit above the critical temperature. To 
perform the simulations at different temperatures we varied [3. To proper 
renormalizc the currents we used nonperturbatively determined .Zy-factors 
for the vector and axial vector channel and tadpole improved Zy-factors 
for the other channels. From spatial CFs we obtained screening masses via 

l/T 

G H {z,T) = j dr J dxdy(J H (T,r)jl(T,r))KCOsh(m%(z-N z /2)). (3) 
o 

They show almost no temperature dependence below T c (fig.2) and the ex- 
pected scaling with the temperature above T c . The screening masses in 
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Figure 3. SF obtained from the MEM at different temperatures in the pseudo-scalar 
channel(left) and in the vector channel(right). For details see text. 



the PS and S channel agree with each other already at 1.5T C which is an 
indication that the Ua{^) symmetry is already effectively restored at this 
temperature. The pole masses in the chiral limit can be addressed by us- 
ing Eq.(l) and a zero temperature ansatz of the SF which leads in the 
pseudo-scalar channel to 

Gps(r T = 0.5, T) f a PS (u;, T)/T 2 T 



T 3 



f ctp S (lj,T)/T* T 
J 1 sh(w/2T) m PS sh(mp S /2T)' [ ' 



This quantity has the advantage to be UV save as opposed to the suscep- 
tibilities xps(T) = / Gh(t,T)cIt. After the extrapolation to the infinite 
volume limit the extrapolation to the chiral limit for 0.6T C and 0.9Tc (see 
fig. 2) leads to a vanishing pion mass which therefore remains a Goldstonc 
boson. The SF obtained from MEM are illustrated in fig. 3 where the SFs 
are only shown up to ui/T = 10 for clarity (see sec. 2). Below T c the quark 
masses are chosen to give an approximately constant pion mass. In this way 
we get a constant vector mass below T c so that we have no indication of a 
change in the rho mass at least up to 0.9T C . Above T c we see the scaling 
of the peak location with the temperature in both channels. 



4. Thermal Dilepton Rates 

In two flavor QCD the differential dilepton rate d t^i is directly connected 
to the vector spectral function av{u>,T) by 

dN a dW 5a 2 (j v (lu,T) 



d 4 xd 4 p dtod^p 277r 2 w 2 (e"/ T -l) 



(5) 
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Figure 4. The vector SF for different cut-offs at 1.5T C and the resulting thermal dilepton 
rate on 16 X 64 3 lattices. 

The Born rate is given by the free vector spectral function <jy ee = 

0. 757. -2 w 2 tanh(cj/4T). To ensure that the results are not influenced by 
lattice artefacts we carried out the MEM analysis on data with different 
cut-offs (sec fig. 4). Employing Eq.(5) to compute the thermal dilepton rate 
we see an enhancement over the Born rate in the region oj /T = 5 — 6 and a 
very close agreement at higher energies 5 . The most striking feature is the 
sharp drop at energies below ui/T = 2 — 3. If this result persist even closer 
to T c this implies that there are no thermal contributions to the dilepton 
rate at low energies which is in contradiction to other approaches 6 . This re- 
sult is also corroborated by a direct inspectation of the correlation function 
at tT = 0.5, Eq.(4), which gives G v (1/2T, T)/T 3 = 2.23 ±0.05/2.21 ±0.05 
at T/T c = 1.5/3. The finitcness of this number already requires a vanishing 
spectral function in the limit oj — > and the fact that it is close to the value 
obtained with the free spectral function, G f v ree (l/2T,T)/T 3 = 2, suggests 
that the contribution from the low ui region is small. 
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